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Two  Exponential  Approximation  Methods 


I.  Introduction 

Two  design  procedures  for  aperiodic,  or  space  tapered,  linear  arrays  are  in¬ 
vestigated  in  this  report  in  a  setting  much  more  general  than  the  usual  setting.  One 
procedure,  due  to  Bruce  and  Unz  [1],  gives  both  element  excitations  (“shadings”) 
and  positions.  The  other  procedure,  due  to  Maffett  [2],  gives  element  positions 
under  the  condition  that  all  excitations  are  unity.  Both  seek  desirable  radiation 
patterns  minimizing  grating  lobes.  These  methods  synthesize  sparse  broadband 
arrays  that  are  less  sensitive  to  frequency  changes  than  periodic  (equispaced)  arrays. 

Using  either  of  these  procedures,  the  designer  must  guess  the  number  of 
elements  required,  perform  the  appropriate  numerical  computations,  examine  the 
resulting  radiation  pattern,  and  then  decide  if  more  elements  are  required  or  if  fewer 
elements  will  suffice.  In  this  report,  error  bounds  are  derived  that  provide  estimates 
on  the  number  of  elements  necessary  for  a  given  degree  of  approximation  of  the 
desired  radiation  pattern.  Thus,  some  of  the  effort  and  guesswork  inherent  in  these 
procedures  can  be  eliminated. 

Neither  of  these  two  methods  is  intrinsically  limited  to  aperiodic  array  design. 
Generalizations  turn  out  to  be  worthwhile  and  of  independent  interest.  Therefore, 
this  report  addresses  only  the  general  setting  from  this  point  on. 

A  complex  valued  function  f  of  a  real  variable  is  defined  to  be  positive 
definite  if  and  only  if,  for  each  integer  n^l,  the  inequality 

n 

Z  f(Xj-X:)  a;  a :  >  0  (1.1) 

i,j  =  l 

holds  for  all  x,,  .  .  .  ,  xn  eR  (the  real  numbers)  and  a,,  .  .  .  ,  an  eC  (the  complex 
numbers).  Bochner’s  Theorem  states  the  following:  If  f  is  a  continuous  function  on 
R,  then  f  is  positive  definite  if,  and  only  if,  there  exists  a  bounded  non-decreasing 
function  V  on  R  such  that  f  is  the  Fourier-Stieltjes  transform  of  V;  that  is, 

rm 

f(x)  =  J  eloxdV(a),  xeR  .  (1.2) 

The  recent  paper  of  Stewart  [3j  gives  references  to  various  proofs  of  Bochner’s 
Theorem  and  its  generalizations.  We  point  out,  for  future  use,  that  (1.2)  im¬ 
mediately  implies  that  the  total  variation  of  V  equals  f(0),  and  for  all  real  x,  f(-x) 
equals  the  complex  conjugate  of  f(x).  Goldberg  [4]  proves  that  any  positive  definite 
function  f  is  such  that  |f(x)|<f(0)  for  all  real  x. 

In  this  report,  we  restrict  our  attention,  for  the  most  part,  to  continuous 
positive  definite  functions  f  on  R  that  can  be  written 

f(x)  =  j  e1®'  dV(a),  xeR,  (1.3) 

J  _  1 
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for  some  real  number  A  such  that  0  <  A  <  °°.  In  other  words,  we  have  assumed  that 
V(cr)  is  constant  for  |o|  >  A.  For  functions  satisfying  (1.3),  we  develop  in  an 
elementary  manner  an  approximation  to  f(x)  of  the  form 

n 

Sn(x)  =  f(0)  2  ^  ak  eiok'  ,  (1.4) 

where  |ak|  <  A,  k  =  I,  .  .  .  ,n.  We  give  an  error  bound  for  this  approximation  in 
Theorem  1.  This  approximation  always  gives  positive  coefficients  and  exponents 
akcR  that  are  located  at  the  roots  of  an  appropriate  orthogonal  polynomial.  We 
suspect  that  these  approximations  are  near-optimal  in  some  well-defined  sense.  (See 
Schabach  (5,  p.  1018)  for  a  relevant  conjecture  about  a  particular  function  f.) 

Under  various  additional  assumptions  concerning  V,  we  develop  an  ap¬ 
proximation  to  f(x)  of  the  form 

Q-w-Ti-”-.  d.» 

where  | erk{  <  A,  k=  1,  .  .  .  n,  and  we  give  an  error  bound  in  Theorem  2.  The  ap¬ 
proximation  Qn(x)  cannot  be  as  efficient  in  general  as  the  approximation  Sn(x); 
however,  Qn(x)  has  the  advantage  of  being  much  more  easily  constructed  in  practice 
for  almost  any  reasonable  n  (say,  n  <  106). 

Note  that  both  the  approximations  Sn(x)  and  Qn(x)  are  readily  written  in  the 
form  (1.3)  and,  therefore,  are  positive  definite.  Hence,  we  must  have 

|f(x)-Sn(x)|  <  | f(x)|  +  |Sn(x)|  <  2f(0),  xrR  ,  (1.6) 

since  Sn(0)  =  f(0).  Similarly,  it  is  always  the  case  that  |f(x)  -  Qn(x)|  <  2  f(0). 

It  will  be  shown  that  Prony’s  method  can  be  used  to  compute  Sn(x).  Although 
Prony’s  method  in  this  problem  must  become  numerically  ill-conditioned  for  n 
sufficiently  large,  it  may  nonetheless  be  useful  for  small  n  (say,  n  <10). 
Numerically  stable  methods  for  computing  Sn(x)  suitable  for  all  n  would  require  an 
algorithm  other  than  Prony’s  method.  This  is  discussed  at  the  end  of  Section  II. 

The  computation  of  approximations  Qn(x)  is  show  n  to  depend  upon  the  ability 
to  compute  the  numerical  value  of  the  inverse  function  of  V  (guaranteed  to  exist  by 
additional  assumptions)  at  specific  points.  The  level  of  difficulty  involved  depends 
on  V,  of  course,  but  the  interval  is  finite,  so  the  problem  seems  to  encounter  no 
inherent  numerical  difficulties. 

An  excellent  bibliography  of  references  to  the  literature  on  exponential  ap¬ 
proximation  is  contained  in  [6]. 

Note  that  if  V  in  (1.3)  is  continuously  differentiable  on  the  interval  (-A,  A),  then 
V'(ff)  >  0,  and 


f(x)  =  J  ela'  V’(o)  da  . 


From  the  Paley-Wiener  Theorem  (see,  e.g.,  [7,  p.  134]),  this  equation  uniquely 
extends  the  domain  of  f  to  all  C  and  that  this  extension  of  f  is  an  entire  function 
of  exponential  type  at  most  A. 


We  close  this  section  with  a  small  collection  of  positive  definite  functions. 
According  to  [3],  Schoenberg  proved  that  fr(x)  =  exp[-|x|r]  is  positive  definite  if 
and  only  if  0  <  r  <  2,  and  Polya  proved  that  any  real,  even,  continuous  function  f 
that  is  convex  on  the  interval  (0,  oo),  that  is,  f((x  +  y)/2)  <  (f(x)  + f(y))/2,  and 
satisfies  limv  f(x)  =  0,  is  positive  definite.  Goldberg  [4,  p.  61]  proves  that  if  f(x) 
is  positive  definite  and  a  >  0,  then  the  function  h(x)  =  f(x)  exp  (-ax-)  is  also  positive 
definite.  Finally,  if  the  function  f  has  a  Fourier  transform  that  is  nonnegative  and 
integrable,  then  the  function  f  is  positive  definite.  Specific  examples  of  functions 
satisfying  this  latter  property  are 


V(+A)  =  1/2A  +  lim0-.A-  V(a),  and  where  I,<x>  is  the  modified  Bessel  function  of 
order  one.  This  function  is  interesting  because,  for  |x|  >  a/A,  it  has  magnitude  not 
exceeding  1,  while  for  |x|  <  a/A,  it  exhibits  very  rapid  growth  achieving  a  maximum 
magnitude  of  cosh(a)  at  x  =  0.  Other  examples  can  be  discerned  in  various  tables  of 
integral  transforms,  such  as  [1 1]. 
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11.  Exponential  Approximation  W  ith 
Arbitrary  Coefficients 


The  idea  developed  in  this  section  for  constructing  approximations  of  the  form 
Sn(x)  is  simply  Gaussian  quadrature.  A  glance  at  equation  (1.3)  reveals  that  we  are 
particularly  interested  in  Gaussian  quadrature  with  respect  to  the  measure  dV(«). 
From  Szego  [12,  p.  25],  a  system  of  orthogonal  polynomials  exists  for  the  measure 
dV(a)  if  V(o)  has  infinitely  many  points  of  increase  in  the  interval  [-A ,  -A]  and  if  the 
moments 

/■+* 

cm  =  J  am  dV(a)  ,  m  =  0,1,2,  .  .  .  (2.1) 

-A 

exist.  Since  V  is  bounded  above,  the  moments  cm  certainly  exist.  If  V  has  finitely 
many  points  of  increase,  then  f  can  be  written  explicitly  as  a  finite  sum  of  ex¬ 
ponentials.  Although  this  special  case  is  not  uninteresting  (in  the  context  of 
economizing  large  finite  exponential  sums),  we  will  avoid  it  by  assuming  that  V  has 
infinitely  many  points  of  increase. 

Let  a,,  .  .  .  ,  an  be  the  abscissas  and  let  b,,  .  .  .  ,  bn  be  the  corresponding  Cotes 
numbers  of  the  n-th  order  Gaussian  quadrature  formula  with  respect  to  the  measure 
dV(a).  Since 

i  bk  =  r  i  dv(o)  =  ao) , 

k  =  I  •'-A 

we  rewrite  the  Cotes  numbers  in  the  form  bk  =  ak  f(0),  k  =  1 . n.  Using  this 

notation,  and  applying  the  quadrature  formula  blindly  to  (1.3)  gives  the  ap¬ 
proximation 

Sn(x)  =  f(0)  2  ake'*k*.  (2.2) 

k  =  I 

where' 

0<ak  <  1,  k  =  1 . n  (2.3) 

a,  +  .  .  .  +  an  =  1  (2.4) 

|«k|  <  A,  k  =  1 . n  .  (2.5) 

These  three  properties  are  immediate  consequences  of  well  known  results  on 
Gaussian  quadrature.  (See  Szeg’d  [12,  pp.  47-49].)  In  addition,  these  properties 
imply  that  Sn(x)  possesses  good  numerical  round-off  error  behavior  w  hen  the  sum  is 
evaluated  numerically. 

We  seek  an  error  bound  such  that 

|f(x)  -  Sn(x)j  <  Rn(x)  ,  xrR  .  (2.6) 

It  is  clear  from  (1.3)  and  the  Riemann-Lebesgue  Lemma  that  f(x)-*0  as  x-*°°.  On 
the  other  hand,  it  is  not  hard  to  see  from  (2.2)  that  Sn(x)  cannot  tend  to  zero  as 

x- *oo.  The  most  that  can  be  expected  is  that  Rn(x)  becomes  "small’’  for  any  fixed  x. 

We  will  show  that  as  n-*00,  Sn  converges  to  f  uniformly  on  any  finite  real  interval. 
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Let  n  >  1 .  For  each  xeR,  let 

£n(x)  =  min  max  |  eio'  -  n,n_,(a)|  *  (2.7) 

-A<»<A 

where  the  minimum  is  taken  over  all  polynomials  n-,n_,(a)  of  degree  at  most  2n-l 
with  complex  coefficients.  We  always  have  £n(x)  <  1  for  all  x,  as  can  be  seen  by 
considering  the  case  5  0  in  (2.7). 

Lemma  1.  Forn^l, 

£„(x)<>/2  —  ,  X£R. 

2-"->(2n)!  (2.8) 

Proof.  From  a  theorem  given  in  [13,  p.78],  for  any  real  valued  function  p(a) 
defined  on  the  interval  [-1,  +1]  and  possessing  n  +  1  continuous  derivatives  on 
(-1,  + 1),  we  have 

Id("  +  !)(£)l 

E  (p)  a  min  max  |  p(a)  -  nn(a)  |  =  - 

-Ko<l  2n(n  +  l)! 

for  some  4,  - 1  <  4  <  +  1 ,  where  the  minimum  is  taken  over  all  real  polynomials  nn  of 
degree  at  most  n.  For  p(cr)  =  cos  crAx  defined  for  a  in  the  interval  [-1,  1]  and  for 
fixed  real  numbers  A  and  x,  we  have 


E;-1  'Pl  =  2^0^  |C°S  £iX| 
^  (Ax)-n  # 
2:n-l(2n)l 

For  q(o)  =  sin  a  Ax  on  [-1 ,  1],  we  have  similarly 

c-  ,  ,  „  (Ax)-n 

E-,n.,(q)<  -  • 

22nl(2n)! 

From  the  definition  of  £n(x),  we  have 


£n(x)  =  min  max|eioi'  -  n,n_,(Aa)| 

-Ko«l 
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<  ieL,(p»  +  eL.khI'1  • 

Substituting  the  estimates  tor  E^p.^p)  and  E2„-i(q)  completes  the  proof. 

We  remark,  but  do  not  prove,  that  an  example  in  Meinardus  [13.  p.  96J  can  be 
extended  and  used  to  show  that  for  fixed  x, 

( A\l-n 

£„(*)  <  _ i _  (1  +  0(1 ) )  ,n— oo  . 

2-"- 1  (2n)! 
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It  seems  reasonable  to  conjecture  that  this  asymptotic  inequality  is  actually  an 
asymptotic  equality.  In  any  event,  we  use  only  (2.8)  in  this  report. 


Theorem  /.  Let  f(x)  be  a  continuous  complex  valued  positive  definite  function  of  a 
real  variable  such  that 


t(x)  =  J  e'“'dV(a)  ,  xeR,  (2.9) 

-A 

where  V  is  a  bounded  non-decreasing  function  having  infinitely  many  points  of 
increase  in  the  finite  closed  interval  [-A,  A],  Then,  for  each  integer  n  >  1 ,  there  exists 

distinct  real  numbers  . on  and  real  numbers  a! . an  satisfying  (2.3), 

(2.4),  and  (2.5)  and  the  additional  condition 


|f(x)  -  Sn(x)|  <  v/T  f(0) 


(Ax):n 
21"-2  (2n)l 


,  xeR,  (2.10) 


where  Sn(x)  is  given  by  (2.2).  Furthermore,  the  left-hand  side  of  (2.10)  is  never 
larger  than  2f(0)  for  all  xeR  and  every  integer  n  >  1 . 

Proof.  Let  a,,  .  .  .  ,  an  be  the  distinct  abscissas  of  an  n  point  Gaussian  quadrature 
formula  with  respect  to  the  measure  dV(a),  and  let  b, . bn  be  the  correspond¬ 

ing  Cotes  numbers.  Let  the  numbers  a,,  .  .  .  ,  an  be  defined  by  the  relationship 
bk  =  f(0)  ak,  k  =  1 ,  .  .  .  ,  n.  Equations  (2.3),  (2.4),  and  (2.5)  are  then  satisfied.  Fix 
xeR.  Let  p*(»)  be  any  polynomial  of  degree  at  most  2n-l  such  that 


max  |ei®'-p*(er)|  =  £n(x)  , 


-A«<KA 


where  En(x)  is  defined  by  (2.7).  Then,  defining  Sn(x)  as  in  (2.2),  we  have 


|f(x)-Sn(x)|  <|f(x)-f(0)I  akp*(ok)| 

k  =  I 

+  |f(0)S  ak  p*(ak)  -  S  (x)| 

k  =  I 


|eio*  -  p*(a)|dV(a) 


ii 

+  f(0)I  ak|p*(ak)  -  eiok'| 


k  =  i 


r‘ *  JL 

<En(x)J  dV(a)  +  f(0)fn(x)Z  ak 

-A  k  =  I 


=  2f(0)£n(x)  . 
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Since  tn(.\)  <  1  is  always  true,  recalling  Lemma  1  completes  the  proof. 


Corollary  1.1.  Any  sequence  of  approximations  S„(x),  n=l,2,  ....  satisfying 
Theorem  1  converges  uniformly  to  f(x)  on  every  finite  interval. 

Proof.  Immediate. 

Corollary  1.2.  If,  in  addition  to  the  requirements  of  Theorem  1 ,  f(x)  is  real  valued, 
then  for  each  integer  n  >  1,  there  exists  distinct  real  numbers  f} 0n  and  real 
numbers  d, . dn  that  satisfy 

0 <  dk  <  1  ,  k  =  I,  .  .  .  ,  n  (2.11) 

d,  +  d2  +  .  .  .  +  dn  =  1  (2.12) 

0  <  /3k  <  A  ,  k  =  1 . n  (2.13) 

and  are  such  that 

|f(x)-f(0)2i  dkcos/>kx|<f«»  (2.14) 

for  all  xeR.  Furthermore,  the  left-hand  side  of  (2.14)  is  bounded  from  above  by 
2f(0)  for  all  xeR  and  every  integer  n  >  1 . 

Proof.  Since  f(x)  is  real  valued,  by  conjugating  (1.1)  we  see  that  it  must  be  even. 
From  f(x)  =  (f(x)  +  f(-x))/2  and  (2.9),  we  get 

f(x)  =  f  cos  ox  dV(ff)  .  (2.15) 

—A 

Furthermore,  the  measure  dV(o)  can  be  taken  to  be  symmetric  about  0.  For  each 
n  >  1,  and  for  each  fixed  xeR,  define 

£j,(.x)  =  min  max  |cos  ax  -  n,n.,(cr)j  , 

-A«a«A 


where  the  minimum  is  taken  over  all  polynomials  Tt-,n_,(cr)  of  degree  at  most  2n-I 
with  real  coefficients.  Hence,  we  always  have  ^n(x)  «  I  by  considering  the  case 
n,n.|(o)  a  0.  From  the  proof  of  Lemma  1 , 


Duplicating  the  proof  of  Theorem  I  with  2n  replacing  n  gives 

2n 

|f(x)  -  f(0)  2  \  cos  okx|  <  2f(0)^n(x) 
k  =  i 

for  the  distinct  real  numbers  o, . o,n  and  real  numbers  a, . a,n  that  are 

the  abscissas  and  Cotes  numbers,  respectively,  of  the  Gaussian  quadrature  of  order 
2n  with  respect  to  the  measure  dV(<*).  These  abscissas  and  Cotes  numbers  satisfy 
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(2.3),  (2.4),  and  (2.5).  Since  the  measure  dV(»)  is  symmetric  about  zero,  it  must  be 
that  a,  =  -a;„,  ~  -a2n-p  etc->  anc*  that  a,  =  a2n,  a:  =  a2n_;,  etc.  Inequality 

(2.14)  follows  immediately  by  taking  dk  =  2alwkand/)k  =  an,k  for  k  =  1, 2 . 

n.  The  properties  (2.1 1),  (2.12),  and  (2.13)  follow  from  (2.3),  (2.4),  and  (2.5).  This 
completes  the  proof. 


Example  1.  The  real  valued  function 

f(x)  =  2  ~^-  =  j'eio'dV(a) 


(2.16) 


with  V(a)  =  a,  -1  <  a  <  1,  is  a  positive  definite  function  on  R.  In  this  case, 
Gaussian  quadrature  with  respect  to  the  measure  dV(a)  is  Gauss-Legendre 
quadrature.  Thus,  from  the  proof  of  Corollary  1 .2,  for  each  n^l.we  have 


sin  x 


0  1 

-Z  dkcos0kx|  < 


*4n 


24n‘2(4n)! 


(2.17) 


where  p ,  Pn  are  the  positive  abscissas  of  a  2n  point  Gauss-Legendre 
quadrature  and  d,,  .  .  .  ,  dn  are  the  corresponding  Cotes  numbers.  This  example 
and  some  computation  provides  one  test  of  the  quality  of  the  error  term.  Let  Rn(x) 
be  the  smaller  of  the  two  numbers  2 f(0)  =  2  and 

x4n  /  xe\4n  f~2 

24n'2(4n).'  \8n/  V  nn  (2,18) 

From  Table  1,  it  appears  that  Rn(x)  is  an  excellent  error  bound  provided 
|xi  «  8n/e.  (Table  1  was  computed  on  a  DEC  VAX  11/780  on  which  the  double 
precision  unit  round-off  error  is  only  4  x  1 0-17.) 

Table  1.  Comparison  of  (2.17)  for  n  =  10  (8n/e  ~  29.43) 


X 

R.oCx) 

f(x)-S10(x) 

maxjf(y)-Sl0(y)|, 

5 

.407  x  10-31 

underflow 

underflow 

10 

.447  x  10-'9 

underflow 

underflow 

15 

.494  x  10-'2 

.491  xlO-'3 

.491  x  10-'3 

20 

.491  x  10-7 

.164  x  10’8 

.  164  x  10"8 

25 

.370  xlO-3 

.290  x  10'5 

.290  x  10-5 

30 

.543 

.664  x  10-3 

.664  x  10-3 

35 

.200  x  10‘ 

.302  x  10-' 

.302  x  10-' 

40 

.200x10' 

.309 

.309 

45 

.200  x  10' 

.501 

.569 

50 

.200x10' 

-.364 

.569 

This  section  is  concluded  by  showing  that  Prony’s  method  (see,  e.g.,  [14,  p. 
378]  or  [15,  p.  340]  )  can  be  used  to  compute  numerically  the  approximations  of 
Theorem  1.  We  need  only  find  the  Gaussian  quadrature  formula  of  order  2n  with 
respect  to  the  measure  dV(a),  which  is  equivalent  to  solving  the  equations 
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cB-  fAan,dV(a)  =  i  bkak"',m  =  0,l . 2n-l  (2.19) 

-A  I'  =  1 

for  b,,  ....  bn  and  Op  .  .  .  ,  an.  Since  ak  must  be  real,  write  sk  =  In  ak,  if  ak  *  0, 
and  sk  =  0  if  ak  =  0.  The  required  equations  can  now  be  written  as 

n 

cm  =  2  bkemsi  ,  m  =  0,1 . 2n-l  . 

k  =  1 

This  form  is  precisely  the  form  required  for  Prony’s  method.  (The  use  of  Prony’s 
method  to  compute  Gaussian  quadrature  formulas  was  pointed  out  to  the  author  by 
Marvin  J.  Goldstein.)  In  principle,  we  require  2n  quadratures,  the  solutions  of  two 
systems  of  linear  equations  each  of  rank  n,  and  the  roots  of  a  polynomial  of  degree 
n  (in  this  case,  all  its  roots  are  known  to  be  real,  distinct,  and  have  multiplicity  one) 
to  compute  one  approximation  for  which  Theorem  1  holds.  Unfortunately,  it  is 
known  [16]  that  any  procedure  that  relies  upon  the  moments  must  become  in¬ 
creasingly  numerically  ill-conditioned  as  n  increases.  Fortunately,  the  use  of 
modified  moments  (i.e.,  replacing  the  am  in  (2.19)  with  some  classical  system  of 
orthogonal  polynomials)  together  with  an  algorithm  other  than  Prony's  method 
often  results  in  a  numerically  well-conditioned  problem  for  finite  intervals.  See  [17] 
and  [18]  for  details. 


III.  Exponential  Approximation  With 
Uniform  Coefficients 


The  idea  developed  in  this  section  for  constructing  approximations  of  the  form 
Qn(x)  in  which  each  exponential  term  enters  the  approximation  with  equal  weight  is 
basically  probabilistic  in  nature.  The  integral  representation  (1.3)  of  f(x)  is  ap¬ 
proximated  by  a  Riemann  sum  whose  subintervals  are  equally  probable  according 
to  the  "probability”  measure  dV(a).  In  this  interpretation,  V(a)  is  a  cumulative 
probability  integral  that  is  used  to  transform  n  uniformly  distributed  points  in  the 
range  of  \'(a)  into  n  abscissas  on  the  real  line  distributed  according  to  the  measure 
dV(a).  (See  [19,  p.  314]  or  [15,  p.  389].) 

Theorem  2.  Let  f(x)  be  a  continuous  complex  valued  positive  definite  function  of  a 
real  variable  such  that 

r‘ l 

f(x)  =  J  el<>sdV(o),  xtR,  (3.1) 

-A 

where  V  is  a  continuous  and  strictly  monotone  increasing  function  throughout  the 
finite  closed  interval  [-A,  A].  Then,  for  each  integer  n  3*  1,  there  exists  distinct  real 
numbers  o,,  .  .  .  ,  an  in  the  open  interval  (—A,  A)  such  that 

|f(x)  -  Qn(x)|  <  2vT  f(0)A|x|/n  ,  xtR  ,  (3.2) 


where 


<3.3, 

Furthermore,  from  the  remark  following  (1.6),  the  left-hand  side  of  (3.2)  is 
bounded  from  above  by  2f(0)  for  all  xeR. 

Proof.  Let  the  real  number  x  be  fixed  throughout  this  proof.  Define,  for  k  = 
0,1,2,  .  .  .  ,  2n, 


uk  =  V(-A)  +  (V(A)  -  V(-A))  k/2n 
vk  =  V-'  (uk)  . 


(3.4) 


Under  the  hypotheses  on  V(o),  it  is  clear  that  V'1  exists  and  is  continuous  and 
strictly  monotone  on  the  closed  interval  [V(-A),  V(A)].  Hence,  the  numbers  vk  in 
(3.4)  are  well  defined  and  are  distinct.  It  will  be  shown  that  inequality  (3.2)  holds  for 


ak  -  v:k-t  ’  k  -  1,2,  .  .  .  ,  n  . 


(3.5) 


Since 


V(vkH)-  V(vk)  =  ( V(A)-V (— A))/2n  =  f(0)/2n  ,k  =  0,1 . 2n-l  ,  (3.6) 

it  follows  from  the  definition  (3.3)  of  Qn(x)  that 
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Q„(M  =  X  f  k  exp(i  v,k_,x)dV(a)  .  (3.7) 

k  =  *  '2k-2 

By  the  Mean  Value  Theorem,  there  exists  £k  in  the  interval  between  a  and  v,k_,  such 
that 

cos  ax  -  cos  v,k  ,x  =  -x(a-v:k _,)  sin  4kx  .  (3.8) 

Thus,  a  <  v-,k_,  implies  a  <k  <  V2k-  ,  and  v->k  ,  <  a  implies  v,k_,  <  £k  <  a.  From  (3.1) 
and  (3.7), 

J-  r':k 

Re(f(.x)  -  Qn(x))  =  2.  J  (cos  ox  -  cos  v,k  ,  x)  dV(o) 

k  =  1  v2k-2 

J2-  rv’k 

=  -X  2.  J  (fl,_v2k-i)sin  £kx  dV(«)  (3.9) 

k  =  1  v2k-2 

and  so,  taking  absolute  values, 

|Re(f(x)  -  Qn(x))|  <  |x|Z  f  ’k  |a-v,k.,|  dV(o) 

k  =  l  v2k-2 
n 

<  |*|Z  (v2k“v2k-2>  (V^v-*k>  ~  V(v,k.,)) 

k  =  l 

=  2A  |x|  f(0)/n  .  (3.10) 

where  (3.6)  was  used  in  the  last  step.  Similarly, 

|Im  (f(x)  -  Qn(x)  )|  <  2A  |.x|  f(0)/n  .  (3.11) 

Clearly,  (3.10)  and  (3.1 1)  together  complete  the  proof. 

Corollary  2.1.  Any  sequence  of  approximations  Q„(x),  n  =  1,2 . satisfying 

Theorem  2  converges  uniformly  to  f(x)  on  every  finite  interval. 

Proof.  Immediate. 

Corollary  2.2.  If,  in  addition  to  the  requirements  of  Theorem  2,  f(x)  is  real  valued, 
then  for  each  integer  n  >  1 ,  there  exist  distinct  real  numbers  p,,  ...  ,  pn  in  the  open 
interval  (0,A)  such  that 

|f(x)-  Z  cosM<«0)A|x|/n.  (3. J 2) 

Furthermore,  the  left  hand  side  of  (3.12)  is  never  larger  than  2f(0)  for  all  xrR  and 
integer  n  >  I . 

Proof.  Recalling  (2.15),  follow  the  proof  of  Theorem  2  with  2n  replacing  n 
throughout.  Half  of  the  resulting  2n  ak’s  are  positive.  Set  the  pf  s  equal  to  the 
positive  ok’s.  The  details  are  immediate.  This  concludes  the  proof. 
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The  proof  of  Theorem  2  requires  that  V  be  continuous  and  strictly  monotone 
increasing.  It  is  not  clear  whether  the  hypotheses  on  V  can  be  weakened.  On  the 
other  hand,  the  convergence  rate  of  the  approximations  Q„  can  apparently  be 
improved  by  making  further  assumptions  concerning  V.  In  general,  however,  better 
than  n _2  convergence  rates  cannot  be  expected.  Consider  Example  1,  where  V(o)=a 
and  f(x)  =  x_l  sin  x.  From  the  construction  indicated  in  Corollary  2.2,  f(x)  is  ap¬ 
proximated  by 


Qn(\)  =  -J-  Z  cos  (2k-!)x/2n  = 


2n  sin  x/2n 


(3.13) 


It  can  be  shown  directly  that 


-QnOOl  < 


lim^n-Ix-'  sin  x  -  Qn(x)]  =  X^n  X 


(3.14) 


(3.15) 


Hence,  in  this  example,  the  correct  convergence  rate  is  precisely  n-2  for  each  fixed  x. 
We  point  out  that  the  upper  bound  in  (3.14)  and  the  limit  (3.15)  follow  by  a  trivial 
application  of  a  suggestive  result  in  Polya-Szego  (20,  Pt.  2,  Ch.  1,  Pr.  11]. 
Regrettably,  their  method  seems  applicable  in  this  application  only  to  the  special 
measure  V(a)  =  a. 

A  further  example  seems  to  indicate  that  the  convergence  rate  of  the  ap¬ 
proximations  Qn  can  lie  between  n_1  and  n~2. 

Example  2.  [4,  p.  22]  Let  A  be  a  finite  positive  real  number.  Define  V  for  non¬ 
negative  arguments  a>0  by 

f  cr(  1  — cr/2A)  ,0<a<A, 

V(a)  =  U/2  ,  A  <  o, 

and  for  negative  arguments  by  the  relation  V(-o)  =  -V(a),  a  >  0.  Thus,  V  is  an  odd 
function  whose  derivative  V'(a)  =  1  -  |a|/\.  Obviously,  for  all  x  #  0, 


f(x)  =  J  cos  a  x  dV(er)  = 


2( 1 -cos  A  x) 
Ax2 


and  f(0)  =  A.  Now  V-1  exists  on  the  interval  [-A,  A]  and,  for  non-negative 
arguments,  is  given  by 

V-'(t)  =  A(1  -  Vl  -  2t/A  )  ,  0  <  t  <  A/2  . 

From  the  construction  indicated  in  Corollary  2.2, 


=  —  Z  cos  l  A(  I  -\/l  -<k-~!:)/n  )  x]. 


(3.16) 
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I  f(x)  -  Q„(x)|  <  A-  |\|/n  ,  xeR  . 


This  estimate  is  not  even  close  to  the  truth.  In  fact,  an  examination  of  Table  2  in¬ 
dicates  that  for  sufficiently  large  n  the  best  error  bound  may  take  the  form 


|f(x)-Qn(x)|  < 


Kx- 
(2n)'  ^ 


(3.17) 


for  some  constant  K.  In  general,  we  speculate  that  if  V  >  satisfies  a  Lipschitz 
condition  of  order  r,  0  <  r  <  I ,  then  the  convergence  rate  is  of  order  I  /n1  +  r. 


Table  2.  Inequality  (3.17)  for  x  =  10,  A  =  1,K  =  102 


n 

Kx2/(2n)2  2 

f(x)  -  Q„(x) 

5 

.316x  10-i 

-.618x10-' 

10 

.1 12  x  10-i 

-.1 17  x  10-' 

20 

.395  x  10-2 

-.153x10-’ 

40 

.140x  10-2 

.649  x  10  4 

80 

.494  x  10--’ 

.163  x  10  4 

160 

.175  x  10-3 

.907  x  lO  4 

320 

.618  x  10  4 

.400  x  lO  4 

640 

.2J8xlO-4 

.160X10-4 

1280 

.772  X  10-5 

.614  X  lO  5 

2560 

.273  x  10-5 

.229x10-5 

5120 

.965  x  IO-« 

.837  x  lO  6 

10240 

.341  x  10-6 

.303  x  10  6 

20480 

.121  x  10-6 

.109x10-6 

40960 

.426  x  lO-1 

.389  x  10-7 

,'r 


IV.  Concluding  Remarks 


The  proof's  in  this  report  depend  heavily  on  the  finite  support  of  the  measure 
dV(a)  even  though  the  construction  of  the  approximations  Sn(x)  and  Qn(x)  can  be 
carried  out  without  modification  on  infinite  intervals  as  well,  provided  V(o)  is 
bounded.  Since  these  proofs  cannot  be  adapted  for  infinite  intervals,  the  ef¬ 
fectiveness  of  the  resulting  approximations  theoretically  remains  an  open  question. 
Intuitively,  however,  it  would  seem  that  only  our  proofs  are  limited  and  that  the 
underlying  approximation  process  is  generally  valid. 

In  computational  practice  the  function  V  is  usually  unknown.  In  many  cases, 
however,  the  given  function  f  does  possess  a  nicely  behaved  Fourier  transform  from 
which  V  can  be  readily  constructed.  The  Fourier  transform  of  f  can,  of  course,  be 
computed  accurately  and  efficiently  in  many  situations  using  fast  Fourier  transform 
(FFT)  methods. 

If  in  the  above,  V  was  not  monotonic,  but  of  bounded  variation  on  R,  ex¬ 
ponential  approximations  can  be  constructed  as  follows.  In  this  case,  there  exist 
monotone  increasing  functions  V  *  and  V  such  that  V  =  V  *  -V  .  For  each  A  >  0, 
define  the  “bandlimiting”  operator  BA  by 

/•* 

BA  f(x)  =  J  ela'dV(o)  .  xrR  . 

Let  1  be  any  finite  interval.  Given  t  >  0,  choose  A  >  0  so  that  ||f-BAf||<£,  where  the 
norm  is  the  uniform  norm  over  the  interval  l.  Now,  let  the  two  functions 

Ba5  f(x)  =  J  e10'  dV  f  (a)  ,  xeR  , 

be  approximated  using  either  of  the  methods  of  this  paper  by  the  two  finite  ex¬ 
ponential  sums,  say,  E  *  (x),  so  that 

llB/f-EMI  <*  • 

Let  E(x)  =  E'(x)-  E  (X)  .  Since  BAf  =  BA‘f  -  BA  f  ,  we  have 

1 1  f  -  E 1 1  =  ||(f-BA  f)  +  (B;f-E‘)-(BA  f-  em!  | 

<  ||f-BA  f||  +  f|  B/f-E*  II  +  II  Ba  f  -  E-||  <  3c  . 

Thus,  exponential  sums  of  degree  not  greater  than  deg(E  * )  +  deg(E  )  may  be 
constructed  to  approximate  f(x)  with  specified  accuracy  on  the  interval  I. 

It  is  well  known  (4,  p.  60]  that  if  fix)  is  measurable,  then  Bochner’s  Theorem 
still  holds  for  almost  all  x.  but  not  necessarily  all  x  as  in  (1.2).  Results  similar  to  the 
results  of  this  report  can  also  be  proven  for  measurable  f  with  careful  attention  to 
certain  details;  however,  this  generalization  is  not  pursued  here.  Similarly, 
Bochner’s  Theorem  has  been  generalized  to  locally  compact  abelian  groups  |4,  p. 
72),  and  perhaps  the  basic  approaches  to  approximation  used  here  can  be  extended 
to  this  much  more  general  setting. 
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We  conclude  by  commenting  that  Bernstein's  Theorem  [21 ,  p.  160)  states  that  a 
necessary  and  sufficient  condition  for  f(x)  to  be  completely  monotonic  on  the  in¬ 
terval  (0,°°)  is  that 

oo 

f(x)  =  f  e°'d  V(a)  ,  0  <  x  <  <»  , 

Jo 

where  V(a)  is  bounded  and  nondecreasing.  It  is  evident  that  the  methods  employed 
in  this  report  can  be  used  in  a  manner  entirely  analogous  to  the  proof  of  Theorem  1 
to  develop  exponential  approximations.  That  is  to  say,  whenever  f(x)  can  be  ex¬ 
pressed  as 

r * 

f(x)  =  J  e°'dV(a) 

Jo 

for  some  finite  A  >  0,  there  exist  approximations  of  the  form 

n 

Tn(x)  =  f(0)X  ake“ks 

k  =  i 

where 

ak>0  ,  k  =  1,  .  .  .  ,  n  (4.1) 


a,  + 


•  •  +  an 


=  1 


(4.2) 


and 


A>  ak >0  ,  k  =  1 . n 


(4.3) 


|f(x)-Tn(x)|  <  2f(0) 


(Ax)-n 

22n 


-Ax/2 


0<  X  <  oo  . 


(4.4) 


An  alternate  approach  for  approximation  of  completely  monotonic  functions  can 
be  found  in  [22]. 
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